Triplet pair correlations in s-wave superfluids as a signature of the FFLO state 
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We show that anti-parallel triplet pairing correlations are generated in superfluids with purely 
s-wave interactions whenever population imbalance enforces anisotropic Fulde-Ferrell (FF) or in- 
homogeneous Larkin-Ovchinikov (LO) states. These triplet correlations appear in the Cooper pair 
wave function, while the triplet part of the gap remains zero. The same set of quasiparticle states 
contributes to the triplet component and to the polarization, thus spatially correlating them. In the 
LO case, this set forms a narrow band of Andreev states centered on the nodes of the s-wave order 
parameter. This picture naturally provides a unifying explanation of previous findings that attrac- 
tive p-wave interaction stabilizes FFLO states. We also study a similar triplet mixing which occurs 
when a balanced two-component system displays FFLO type oscillations due to a spin-dependent 
optical lattice. We discuss how this triplet component can be measured in systems of ultra-cold 
atoms using a rapid ramp across a p-wave Feshbach resonance. This should provide a smoking gun 
signature of FFLO states. 



PACS numbers: 67.85.-d,03.75.Ss,71.10.Pm,74.45.+c 

Conventional (homogeneous and unpolarizcd) BCS 
states survive in the presence of a Zeeman field as long 
as the condensation energy is larger than polarization 
energy gain in the normal phase |Tj. A deep argumenta- 
tion of why these two are the only allowed homogenous 
and isotropic configurations for usual s-wave supercon- 
ductors can be found in Ref. [2]. Fulde-Ferrel (FF) [3] 
and Larkin-Ovchinnikov (LO) [4] found that a compro- 
mise between BCS condensation and polarization could 
occur in a region of the phase space lying between the 
unpolarizcd BCS and the polarized normal Fermi gas 
phases. Those states are characterized by the Cooper 
pair wave function having one (FF) or two (LO) Fourier 
components of nonzero total momentum. As a result, ro- 
tational (FF and LO) and translational (LO) symmetries 
are spontaneously broken. 

The FFLO state is expected to occur in a variety of 
contexts including heavy fermion systems 5 , high den- 
sity QCD [5] and neutron stars [7J. However, there is 
much controversy on the existence of FFLO states in su- 
perconducting systems [8]. Thus, there is a strong mo- 
tivation to explore this intriguing state of matter in sys- 
tems that allow detailed control of their parameters, such 
as ultracold fermions [3]. Moreover, cold atomic gases 
can be used to implement experimental probes that are 
very different from what is available in electron systems. 
Hence, experiments with ultracold atoms can provide 
unique signatures of exotic many-body states. The goal 
of this paper is to show that an unequivocal signature 
of the FFLO phase can be obtained through direct mea- 
surements of the triplet pairing component in the FFLO 
phase in cold atomic gases. In this paper we show that 
triplet correlations always accompany FFLO-type states. 
These correlations are very difficult to probe in electron 
systems (see however |10j and references therein). On 



the other hand, we show how a variation of the tech- 
nique of rapidly ramping along a (in our case, p-wave) 
Feshbach resonance gives a direct access to these triplet 
correlations. This method can be used for an indepen- 
dent test of FFLO physics beyond those already proposed 
[18H21) . Our method can be implemented with current 
experimental techniques [2"!H 123j . 

Another motivation for our analysis comes from recent 
theoretical studies, which found that p-wave attractive 
interactions are surprisingly efficient in stabilizing the 
FFLO state. Even when the p-wave interaction is but 
a small fraction of the dominant s-wave attraction, it 
considerably enhances triplet pair correlations at temper- 
atures well above the nominal p-wave attraction critical 
temperature. Such behavior appears in low-dimensional 
superconductors under a parallel magnetic field (which 
yields a Zeeman but not an orbital effect) [ITlfTS] and 
dipolar cold gases [14]. Triplet correlations in the pres- 
ence of attractive p-wave interactions seem to reside in a 
narrow band of Andreev states, which also host the po- 
larization and which are localized near the nodes of the 
s-wave gap [TH1 [H] ■ 

Earlier analysis [ITHl"4"l [T6] was based on the assump- 
tion that the presence of an attractive p-wave channel 
is necessary to find triplet pair correlations (see however 
[T7] for the case of Abrikosov lattices). In the present 
work we argue that a discussion centered on the gap 
tends to miss an important property of the Cooper pairs, 
namely, that triplet correlations are independently pro- 
duced, even without the help of interactions in the p-wave 
channel. In fact, we show that an inhomogeneous s-wave 
order parameter together with polarization is sufficient 
to generate a triplet component in the Cooper pair. The 
effect of any residual p-wave interaction is to lower the 
energy of the FFLO state using the already macroscopic 
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occupation of the triplet component. 

We consider two ways of achieving FFLO type states, 
spin imbalance and spin-dependent optical lattices [25] . 
For simplicity of presentation we focus on ID scenar- 
ios of neutral two-species Fermi superfluids that show 
a competition between traditional s-wave pairing and 
paramagnetic order without the intervention of the mag- 
netic orbital (Meissner) effect. Importantly, the inter- 
action between species reduces to an s-wave (singlet) 
channel, namely, guy / dx^(x)^(x)^ i(x)^>^-(x) , where 
<?id is the effective ID coupling constant. We will use 
ginkp/EF — —2.04 as in Refs. [25l[26], corresponding to 
an interaction to kinetic energy ratio of — m<?iD /h 2 = 1.6. 

It has been recently noted [2TH5U] that, in the nor- 
mal (non-BCS condensed) phase, the intra-species inter- 
action can be strongly modified (and in particular turned 
attractive) by the polarization of the medium. As a con- 
sequence, at sufficiently low temperatures, normal phases 
are unstable against the formation of intra-species p-wave 
superfluids. In the unpolarized case considered in Ref. 
|28j . the s-wave interaction always dominates over the 
effectively induced parallel-spin p-wave attraction. We 
conjecture that in such a case, once the s-wave superfluid 
is formed, the resulting reduction of the spin fluctuations 
(due to the formed Cooper pairs) will lower the effective 
parallel spin p-wave attraction more than envisaged in 
Ref. [35]. Thus we have chosen to work with systems 
where the induced p-wave attraction can be neglected. 

To develop a microscopic description of the system, 
we introduce the Bogoliubov transformation ^ a {x,t) — 
Y Jk \ u ^)e- i ^ t c ka +Gvl^x)e l ^ t cl s \ where a := -a, 
c and denote the quasiparticles and the sum should be 
restricted to avoid double counting (we use uj k(J > with 
k the quasiparticle index composed of a quasimomentum 
and the band index), which yields the Bogoliubov - de 
Gennes (BdG) equations for the quasiparticle wave func- 
tions: 
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H a A(x) ' 
A*(x) -H s 

h 2 d 2 

2m dx 2 



Uka(x) 
Vk<j{x) 



Uka(x) 
Vka{x) 



+ V„{x) - +gir>ria(x), (1) 



where we have introduced the one-species density, 
n a (x) := n a (x,x), with n a (x,x r ) := ( i $i T (x)^ cr (x')) the 
one-particle density matrix (normal average), and the 
s-wave gap, A(x) := —gmF(x,x), with F(x,x') :— 
( i Sf^(x)'^/-f(x')) the pair wave function (anomalous aver- 
age) [3 1J . These normal and anomalous averages read 

n t (x,x') = ^f{uj k ^)ul^{x)u k ^{x') 
H(x,x') = ^2f(-Ukt)v k t{x)v^(x') 
F(x,x') = -J2f("k\K t (x)u kt (x% (2) 
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where f{oJ ka ) = [1 + exp(/iaj/ co -/fcsT)] _1 , and are to be 
solved self-consistently. Notice that these expressions 
only contain u krT ,v ka ,uj ka - for a =f an( l t ne sums are 
unrestricted [2"6] . 

It can be shown (see Ref. [2]) that the only 
macroscopically occupied eigenstate 4> a p(xi, ^2) of the 
two-particle density matrix is trivially related to 
the pair wave function via a scaling F(xi,x 2 ) = 
\ZNo<Plt(xi,x 2 ), where <f> a& (x,x') = -<j) aa {x' ,x) is 
the wave function of the Cooper pairs, normalized as 
/ dx\d,X2\4>a,p(%ii x 2)\ 2 = 1, and Nq is the number 
of Cooper pairs. Its singlet/triplet components are given 
hy F s / t {x 1 ,x 2 ) = -[F(x u x 2 )±F{x 2 ,x 1 )]/2. Shifting to 
the Cooper pair center-of-mass and relative position co- 
ordinates, X := {x\ + x 2 )/2, y := x\ — x 2 , F s ' t (X; y) := 
F s l l {X + y/2,X — y/2), and using fermionic permuta- 
tion symmetry, it is immediate to see that F t (X; 0) = 
and dF s (X;y)/dy\ y=0 = 0. The quantity F s (X;0) = 
A(X)/<7id is directly related to the conversion into s- 
wave molecules when ramping through an s-wave Fesh- 
bach resonance (see [25]). In this paper we will demon- 
strate that G t (X) := dF t (X;y)/dy\ y= o measures the 
conversion to p-wave molecules when ramping through 
p-wave Feshbach resonances (see below). 

We start by considering spin imbalanced systems with- 
out spin-dependent potentials, V rj {x) — 0. A simple ex- 
planation of the emergence of triplet pair correlations can 
be attained by studying Eq. ([I]), for a =f, in the partic- 
ular case where the Hartree terms are ignored, with the 
chemical potential difference playing the role of an effec- 
tive exchange (Zeeman) field h \= — As a result, 
H a ->■ H + ah/2, where H = -(h 2 /2m)(d 2 /dx 2 ) - fi 
with fi = (/i-j- + ^i)/2. The resulting eigenvalues are 
^fct = £ k + h/2, with e k the eigenvalue for zero exchange 
field. This problem shows a further property which is 
common to unpolarized systems [32] : given a e k solu- 
tion, Xk(x) '■= ((u k (x), Vk(x)) T , there is another one 
with eigenvalue — e k , and wave function —i<J y x k (x) — 
(-v* k (x),u%(x)) . The identity can then be written: 

S(x-x')= (Xk(x)xi(x') +(Tyxl{x)xl{x')(T l 

e k >Q 

(3) 

After some algebra, we find for the anomalous sin- 
glet /triplet pairing components and polarization at zero- 
temperature: 

F s {x,x') = -- [v* k (x')u k (x)+v* k (x)u k (x')} 

e k >h/2 

F\x,x') = -~ J2 K( x >k(x) ~ vt{x)u k {x')\ 

Q<<L k <h/2 

p(x) = - l\u k (x)\ 2 + \v k (x)\ 2 }. (4) 

0<e k <h/2 

where p(x) = n^(x) — n^(x). Equation Q clearly reveals 
that triplet correlations and spin polarization are closely 
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related because they are supported by the same set of 
quasiparticle states. It is known [15] [33] that Andreev 
sub-gap states exist (in the LO case) around the nodes 
of the gap. According to the present argument, they 
contribute to polarization and triplet correlations only 
when a Zeeman field is present. 

A fully self-consistent computation of a LO state is 
shown in Fig. [TJ The triplet-condensate fraction for this 
case is / dXdy\F t (X; y)\ 2 /N = 0.42. A clear correlation 
is observed between the maxima of the polarization p(x) 
and those of the triplet correlation profile G*(x), both 
peaking near the nodes of the s-wave gap. 
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2X/X 

FIG. 1: The polarization (n-f(X) — n\.(X))/kF in blue dashed- 
dotted, pairing A(X)/Ef in red and triplet order parameter 
G t (X)/k F in green, as functions of 2X/X for a LO state with 
a global polarization (n^ — 714.) = 0.162&F /tt, &fA = 20, T = 
and gmkF/EF = —2.04. 

We now discuss a different route to FFLO type states 
based on a spin-dependent optical potential [25] V\(x) — 
— V\,(x) = Vq cos(27ra/A). Cooper pairs experience 
stretching forces which peak at the zeros of the poten- 
tial, i.e. each spin species is attracted to its respective 
minimum. The resulting elongation of the pairs gener- 
ates local triplet p-wave mixing and a lowering of the 
gap amplitude. As Vo increases, 7r-phases (zeros of the 
gap) appear at some critical values of the potential ampli- 
tude. Then the triplet weight is displaced and centered 
around the new gap nodes, changing the configuration 
abruptly. This results in a discontinuous behavior of the 
total triplet and singlet pair components, as can be seen 
in Fig. [2j Successive jumps take place when a transition 
to a higher order 7r-phase occurs. Clearly, the zeros of 
the gap favor the appearance of triplet correlations. 

In Fig. [3] we show the absolute value of the Fourier 
components of the triplet order parameter. These quan- 
tities are important because, as we show below, they di- 
rectly correlate with the number of p-wave molecules ob- 
tained after a ramp. Specifically, to observe the triplet 
correlations predicted in the present work, we propose a 
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FIG. 2: Singlet (/ dydX\F s (X;y)\ 2 /N , blue) vs triplet 
(/ dydX\F t (X;y)\ 2 /No, red) condensate fractions for a ir- 
phases structure in a spin-dependent lattice potential of wave- 
length A = 30 as a function of the optical lattice strength 
Vo/Ef, for interaction strength giukF / Ef = —2.04 and zero 
temperature. 



combination of a 7r/2 RF-pulse followed by a ramp over 
a p-wave Feshbach resonance, which is a variation of a 
well known technique [3"5Tf3T)] . 
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FIG. 3: 7r-phases setup plot showing the absolute value of the 
Fourier components, G m '■= G <3x= 2 7 rm/A from Eq. dfjl, of the 
triplet order parameter, as function of the lattice potential 
strength Vo (parameters as in Fig. [2|. Only odd components 
are non-zero for Vo/Ef < 0.39, whereas only even ones are 
non-zero for 0.39 < Vo/Ef < 0.69, and so forth (opposite 
behavior as the gap, see [25|h 

In general, p-wave Feshbach resonances can be found 
both between identical or between different hyperfine 
species. However, in the case of 40 K [30] there is no 
interspecies p-wave Feshbach resonances yet found. 6 Li 
does have those resonances, but very far away from the 
834 Gauss s-wave resonance (or from the other s-wave 
Feshbach resonances). However 40 K shows a double p- 
wave Feshbach resonance among bb hyperfine states at 
198.4 and 198.8 Gauss, very close to the 202.1 Gauss s- 
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wave Feshbach resonance. We might use a tt/2 RF-pulse 
to partially convert some of the Cooper pairs (made of 
different hyperfine species) into homo-species pairs and, 
before any relaxation takes place, ramp over a p-wave 
Feshbach resonance (between same hyperfine species) to 
convert them into p-wave molecules. 

To estimate the number of p-wave molecules (of bound 
energy e = —h 2 /map and center of mass momentum Hq) 
formed it is enough to compute the average of the number 
of molecules operator in the state before the ramp (fol- 
lowed by the 7r/2-pulse, Ut/^fU^ = (tf t + <^)/V5, if 
needed) [551155]: 

n-qe '■— y ] (b qem bqem) 
ro=-l,0,l 

r e iq-(r 1 +r 2 )/2 

b lem ■= 2^ d 3 ri d 3 r 2 -j= e m(ri-r 2 ) Q £ 

*l,(ri)*J(r a ), (5) 

here <p e m{r)ap denotes the spinor a, /3 components of the 
molecular orbital p-wave function of angular momentum 
projection quantum number m — —1,0,1 and energy 
e and L 3 is the volume. Assuming quasi-one dimen- 
sional configuration ^(ir) = J d 2 p(f){p)^\(x : p), where 
4>{p) — exp[— p 2 /(2a 2 L )]/(^/Tra±) is the transverse con- 
finement wave function and p refers to the transverse 
coordinates, r — (x,p), the averages can be factorized 
as (*4,(ri)* t (r 2 )) = 4>{p 1 )4>{p 2 )F{x 1 , x 2 ), the same oc- 
curring to the one-particle density matrix terms. Surviv- 
ing terms from the one-particle density matrix n a (xi, x 2 ) 
can be shown to lead to a structureless continuous back- 
ground. The reason is that they always contain one co- 
ordinate from each molecular state involved. The other 
terms, which come from the pair wave function before the 
RF pulse, conspire to depend on its triplet part alone and 
each depends on one set of molecular state variables (s- 
wave components of the Cooper pair are lost due to con- 
traction with the p-wave symmetry molecular state). The 
length scales of variation of F(X; y) as a function of y are 
kp and the BCS coherence length £ = Hvp/Ao, with Ao 
the zero temperature gap. We assume that the the other 
scales are much shorter, namely kpa±, kpa p <C l,kp^. 
Then, we can Taylor expand F t (X;y) in the y variable 
and take the lowest non-zero term, which by symmetry 
is the gradient F t (X; y) ~ yG t (X). The final result is: 

n qx = L\G q J 2 16'ja p a 2 L 

G qi := I J dXe i ^ x G t {X), (6) 

where 7 is a dimensionless factor which can be further 
estimated as 7 ~ 3 for a p a±, and 7 ~ (a p /a±) 4 /3 for 
a p <C aj_ ■ In the case in which the p-wave Feshbach res- 
onance allows for interspecies molecular formation, there 
is no need of tt/2 pulse and the results are the same but 
with a signal four times bigger. 



In the case [UJ of 40 K, with a ID geometry of tubes 
that are optically trapped with a superposed magnetic 
field to control the interaction via the Feshbach res- 
onance at B = 202. 1G and N ~ 100 per tube of 
length L ~ 40/im, we have n ~ 2.5 pm~ 1 and kp = 
Txn/2 ~ 3.93 /iin -1 . For simplicity wc neglect the ex- 
ternal confinement. Using a± ~ 60.3nm, a p ~ 132. 3nm 
and g 1D kF/EF = -2.04 wc find n^^/A ~ 0.09 for 
G fe ~ 0.02fc|, with momentum q x ~ 1.23/im, like in Fig. 
[T] However, the number of tubes in [IT] was 4900 so 
we expect a signal of 440 molecules, which should be de- 
tectable. 

The effect of induced parallel-spin p-wave interactions 
[27H50] on the symmetry of the Cooper pair wave func- 
tion can be neglected at temperatures higher than the 
transition temperature to the p-wave superfluid. On 
the other hand, a non-FFLO state with sufficient p- 
wave character would not show two peaks (with opposite 
center-of-mass momenta) in the measured distribution of 
detected molecules. Rather, it would show a single peak 
with zero center-of-mass momentum. Thus the proposed 
detection scheme can be used to rule out such non-FFLO 
states. 

For the sake of simplicity we have focused on one- 
dimensional situations. The analytical results presented 
here are however not restricted to ID systems. In partic- 
ular, the generalization to higher dimensions of the con- 
nection between polarization and triplet-pair correlations 
(see Eq. Q is straightforward, although the geometrical 
configuration of the s-wave gap profile can be consider- 
ably more complex. 

In summary, we have shown that there exists a very 
general and close relation between triplet component and 
polarization, which explains the findings of [161 117j . Our 
analysis also shows why even a small interaction in the 
p-wave channel (with antiparallel spins) can have strong 
effect on the the stability of the FFLO phase. When 
the FFLO state is formed, its energy will be lowered by 
the attractive p-wave interaction component. This triplet 
component condensation energy provides a common and 
simple explanation of the phenomena mentioned in |11| - 
03]. Similarly, a repulsive p-wave interaction should re- 
duce the stability of the FFLO phase, as suggested in 
some of the aforecited works. 
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